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Abstract 

■ The abelian projection in SU(N) gauge theories is discussed in detail, as well as the 

construction of a disorder parameter to study dual superconductivity as a mecha- 
nism for color confinement. If the ideas of the large N limit are correct, a universal 
iV-independent behavior is expected for the suitable rescaled disorder parameter as 

C*~) . a function of A = g 2 N. 
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1 Introduction 



There exists evidence from lattice simulations [1-3] that color confinement 



D 
> 

X 

is produced by condensation of magnetic monopoles in the vacuum, i.e. dual 
superconductivity. Vacuum behaves as a dual superconductor in the confining 
phase, and goes to normal at the deconfining phase transition. 

The evidence refers to SU{2) and SU(3) pure gauge theories. Preliminary 
data indicate that the same mechanism is at work in QCD with dynami- 
cal quarks [4], a fact which is in line with the ideas of iV c — >• oo [5,6]. As 
N c —* oo with A = g 2 N c fixed the theory should preserve its structure: correc- 
tions 0(1/N C ) are expected to be small and under control. If this is true the 
mechanism of confinement should be the same at all values of N c and also in 
full QCD, quark loops being non-leading in the 1/N C expansion. 
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A direct check of these ideas can be done by exploring, by the same tech- 
niques used in Ref. [1,3], the symmetry of the confining vacuum in SU(N) 
theories. The technique used there was to measure the vacuum expectation 
value (fj) of an operator //, which creates a monopole as a function of the 
temperature, across the deconfining phase transition. A non zero (//) implies 
dual superconductivity. This is exactly what is found below the deconfining 
temperature T c . Above T c , (fi) vanishes, and sectors with different magnetic 
charge are superselected. 

Monopoles are defined by a procedure called abelian projection [7], which we 
summarize in the case of the SU (2) gauge theory. 

Let = • a/2 be any field in the adjoint representation, and = 0/|0| its 
direction in color space, which is defined everywhere in a configuration except 
at zeros of 0. 

A gauge invariant, color singlet field strength tensor can be defined: 

Ffj, u = • G>„ - • (D^ A D u 4>) (1) 

where 

= d^A u - d u A^ + gA^ A A v 

is the color field strength tensor, and 

= - gA^)i 

is the covariant derivative of 0. The two terms in Eq. (1) are separately color 
singlets and gauge invariant. The combination is chosen in such a way that 
bilinear terms in A fl A u , and A^dycj) cancel. Indeed, by explicit computation: 

F\iv = <9 M (0 • A v ) - d v {4> ■ A,) - • (dj A dj) (2) 

If we gauge transform to make = const, e.g. = (0, 0, 1) then the second 
term in Eq. (2) vanishes and becomes abelian 

F, u = d,Al - d v A\ (3) 

This gauge transformation is known as abelian projection on 0, and is defined 
up to a residual U(l), corresponding to rotations around 0. 
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Monopoles can be present in a noncompact formulation of the theory when 
d^F* u can be different from zero [8]. A non zero magnetic current exists 

3u = d»F; v (4) 

which is identically conserved 

Vj v = (5) 



The corresponding U(l) symmetry can be either realized a la Wigner or Higgs 
broken. In the first case the Hilbert space is made of superselected sectors with 
definite magnetic charges. In the second case, under very general assumptions, 
the vacuum behaves as a dual superconductor. 

The detection of dual superconductivity has been successfully performed for 
SU(2) and SU(3) gauge theories [1-3]. In this paper we want to analyze the 
abelian projection and the construction of the disorder parameter \i for generic 
SU(N). 



2 Abelian projection in SU(N) gauge theory 



In analogy to the construction for £77(2), for an arbitrary operator <fi in the 
adjoint representation 



= E <t> a T a 



we can define a field strength tensor F^. Here T a (a = 1 . . . N 2 — 1) are the 
generators of SU (N) in the fundamental representation, with normalization 



Tr 



rparpb 



-5 ab 
2 



We write 



F^ = Tr{0G^} - ^Tr{0[/J M 0,/J,0]} (6) 



The normalization of 4> has been left indeterminate. A change in the nor- 
malization reflects in a change of the relative coefficients of the two terms in 
Eq. (6). The notation is the usual one 
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= A«T a 

G>„ = dfj,A v - d u A /I + ig [Ap, A v \ 
£> M = <9 M -ig [A^<p] 

We want to investigate for what choice of 0, if any, bilinear terms in A^A U , 
A^d^cj) cancel in Eq. (6). If this happens 



F IM , = cl Tr {<f>{d»A v - cU M )} - -Tr {0[<9 M 0, d v <t>]} 



(7) 



and Fpv becomes abelian in the gauge where is diagonal, 

F iiV = cYx{d^A v )-dMA { ,)} 



(8) 



The condition for the cancellation of bilinear terms A^A U in Eq. (6), for any 
A I1 A U) reads 



X 2 (A„ 0) = Tr (0L4 M , A v \) + Tr (0 [[A„ 0], [A v , 0]]) = 



(9) 



If 0o is a solution of Eq. (9), then 

<j ) = U(x)<j )0 U^x) 



(10) 



with arbitrary unitary matrix U (x) is also a solution. Indeed 

X 2 {A ti {x),4>) = X 2 {U\x)A lx {x)U{x),^) = 



since Eq. (9) holds for any choice of A^. In particular one can choose 0o diag- 
onal. The generic 0o, diagonal and constant, obeying Eq. (9) is (see appendix 
A) 



^0 



diag 



p_ p_ p_ _q_ _o_ _o_ 

N'N'"''N' N' N''"' N 



V 



11^ 



where p + q = N , and q — 1 . . . N — 1. If the solution is continuous, it must 
be constant, since it cannot jump from one of the solutions in Eq. (11) to 
another. So the generic solution to our initial problem, which we call $, is 
given by Eq. (10), where O can be any of the matrices $q °f Eq. (11). 

Let us now show that, if X 2 = 0, then the terms A^d u & also vanish. Let us 
call Xi(An, $) such terms, we have 
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Xi(A, $) = Trmd„A u - d u A^)} + Tr{$ [[A„, $], d v $]} 

+Tr{$ [<9 M $, [A,, $]]} - Tr{d^A v ) - d^A^)} 
= Tr{$ [[Ap, + Tr{$ [d^, [A,, $]]} 

-Tr{c^A, - d„cM M } (12) 

Now [9] 

and hence from Eq. (12) 

Xi(A, $) = Tr {$ $], [Q u , $]] + $ [[fy, $], [A,, $]] 

(13) 

which can also be written 

x t (A, $) = x 2 (a + $) - x 2 (n, $) - x 2 (a, $) 

Since Eq. (9) is valid for arbitrary A^, one has: 

X 1 (A,$)=Q if X 2 (A,$)=0 



Then 



Ffiv = Tr{d^A v ) - d^AJ} - l -Tr{$[d^,d u <S>}} (14) 



In the abelian projected gauge, where $ = $q, 

F, v = Tr R(<t>A„) - (15) 

On the other hand, the 't Hooft tensor for a pure gauge field uj^ = —i (d^U) W 
is zero 

= F^u) = Tr{«9 M ($^) - «9,($^)} - ^Tr{$[9 M $,^]} (16) 

and finally, by use of Eq. (16), Eq. (14) can be rewritten 

F, u = Tr{d^A v ) - d^AJ} - Ti{d^u u ) - d^u,)} (17) 

showing that F^ u obeys Bianchi identities. 
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In conclusion an "abelian" field strength F q can be defined for each of the 
fields <& q of Eq. (11), or any other field related to them by a gauge transfor- 
mation. The normalization of the solution is fixed. The invariance group of $o 
is 

SU(q) x SU(N -q)x U(l) (18) 

and Q q (x) belongs to the coset of $q 

^ q (x) = U{x)$lU\x) (19) 

$o defines a symmetric subspace [10], in the sense that the full algebra of 
SU(N) is the sum of the subalgebra L of the little group of $q pl us the 
complement U to the full algebra, and 

[L , L ] E L [L', L ] e L' [I/, L'\ E L (20) 

The main property of symmetric spaces is that any element of the group U 
can be uniquely split in the form 

U = e iL 'e iLo (21) 
a property which will be used below. 

Let us now discuss the abelian projection in the light of the above results. 
Following Ref. [7], let us consider a generic operator X, that transforms co- 
variantly under gauge transformations. X can be diagonalized by a gauge 
transformation: 

X D (x) = U x {x)X{x)U] c {x) (22) 

For each $q in Eq. (H), a field transforming in the adjoint representation of 
SU (N) can be defined: 

&(x) = rf x (x)& U x (x) (23) 

These fields $ 9 can now be used to define N — 1 gauge-invariant field strength 
tensors: 

F* = Tr{$^} - -Tr{^[d^,d v ^]} (24) 
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In the gauge where X is diagonal, & q (x) = $q and, according to the above 
results, these tensors reduce to the abelian form: 

F% = d,Tr {& A„} - d u Ti {<&qA^} (25) 
The diagonal matrices <3>o form a complete set of diagonal matrices, and hence 

N-l 

X D (x) = c ^) $ o ( 26 ) 

q=l 

We denote by a a the diagonal matrices associated to the simple roots via: 

a a = a* H h 

If the generators of the Cartan subalgebra Hi are written in the standard 
basis: 

/ m \ ^ 

{H m )ij = ^2 SikSjk — m Si tm+ i5j tm+ i -. 

\k=i I yj2m(m + 1) 

or 

m 

# m = diag(l,l,i...l, -m,0...0)- ' 



^2m{m + V 

the matrices associated to the simple roots a a are: 



a 1 = a^ = ^diag(l,-l,0,...,0) (27) 
a 2 = a^ = idiag(0,l,-l,0,...,0) 



a,a+l 

a a = a°Hi = - diag (0, 0, ... , 0, l^T, 0, . . . , 0) 



It is then trivial to check that 



Tr {^ a b } = U 



and therefore 



ab 

2 



Tr {X D (x)a a } = ^c a (x) = X a D - X a D +l 



f28) 
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For each point x such that c a (x) = 0, two eigenvalues of X become degenerate: 

X%{x) - X a D +1 = 0, 

and the gauge transformation Ux, Eq. (22), becomes singular. Such a singu- 
larity behaves as a magnetic charge with respect to the U(l) group of eq.(18) 
for q = a. 

Let A® be the diagonal part of the gauge field in the abelian projected gauge [7] 

<=diag («;,...,<) (29) 

The diagonal matrix can be expanded in the form: 

K = ^> % (30) 

The N — 1 abelian photons coincide with the abelian fields defined via the 
abelian projection: 

^{KK]=K ( 31 ) 

The singularities of the gauge transformation Ux at c a (x) = are magnetic 
charges with respect to the abelian field a" a result that will be used to 
construct magnetically charged operators in the next Section. 



3 Construction of the disorder parameter. 

In the abelian projected representation, in which the operator $ is diagonal, 
the generic link U^{x) can be cast in the form 

U,{x) = V^\x)Df{x) i = \...N-\ (32) 

where D % = exp(ic^o;*), a 1 being the diagonal matrix corresponding to the 
simple root a 1 . 

Indeed since the algebra of the little group of <fi a , L , defines a symmetric space, 
Up can be uniquely split as U^ = e lL 'e tL °, as in Eq. (21). a 1 is an element of 
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L and commutes with all the others, since a % is invariant under e tL ° 

Tr {$ a a 6 } = ^5 ab = Tr [e iL » <$> a e~ iL ° a b ) = Tr [<5> a e~ iU) a b e iL °) (33) 

a 1 is the generator of the subgroup U(l) of Eq. (18). Therefore 

= e^'e^e^' (34) 

where L is L minus a 1 . The plaquette 11^ (x) = U^{x)U v {x + p)UUx + 
v)Ul{x) can then be rewritten as the product of matrices of the form in 
Eq. (32) 

IV(:r) = V,(x)D^x)V u (x + /}) + £)Z>t (x + ^(x + £)Z>t (x)l/t(x) 
= ^(x)K(x + /^(x + ^(x^x) A,(* + fi)Dl{x + f>)£>i(x) 

= n^(x)n^(x) (35) 

An abelian U(l) plaquette is thus defined. However an alternative way of defin- 
ing the abelian plaquette would be to operate the separation of U^ v directly 
as done for the single links. It is easy to see that the two definitions differ by 
terms 0(a 2 ). Indeed the second definition can be obtained by factorizing e lL ° 
from the product V l V in Eq. (35). The resulting Lq would come from higher 
terms in the Baker Haussdorf formula, which are 0(a 2 ) and higher. 

The lattice abelian projection is therefore intrinsically undefined by terms of 
order 0(a 2 ). A similar ambiguity comes out if the abelian field is defined by 
2 x 2 or 2 x 1 Wilson loops, instead of the plaquette. 

U(l) monopoles are defined as in Ref. [11]. The angle 9^ is defined by the 
equation 

U% = (36) 

The magnetic current defined as 



3» = K0% 
1 
2 



jfj, identically vanishes by the Bianchi identity. However, if the angle 9^ v is 
defined modulo 2n: 



■■ Ofu, + 2-jrn^, -7r < 9^ < 7r (38) 
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then the magnetic current 



U = (39) 

can be different from zero and is conserved, A^j M = 0. The term propor- 
tional to counts the Dirac strings going through the plaquette, which are 
invisible. 

A monopole at a fixed time exists in an elementary spatial cube such that one 
of the faces has = 1, the others = 0. The visible Q* has then a flux of 
2n, which is balanced by the outgoing invisible string. 

The operator ji which adds a monopole at the site y and time y° to a generic 
configuration, can be constructed as follows 

li = exp {-(3^(Soi(yo) ~ S 0i (y ))} (40) 



Soi are the terms of the action involving space-time (0,i) loops, with a space 
link at yo and the others at t > y . For example for the Wilson action 

S 0i = J2 u oi(yo,n) 

We will recall the construction for Wilson action: the generalization to actions 
containing loops other than the plaquette, e.g. improved actions, is straight- 
forward. 

In the favored abelian projection, according to Eq. (34) 

n i(yo, ft) = n i(z/o, n) n w(s/o, ft) (41) 

S is obtained from S by the following substitution in U 0i (y ,n): 

U t (y , ft) - Uiiyo, fi)e iaab ^ = U'M, ft) (42) 



where bf- is the transverse component of the vector potential generating at 
n a monopole sitting at y, dibf = 0. Classical gauge ambiguities in bt are 
contained in the longitudinal part, and do not affect the definition Eq. (42). 

The plaquette U 0i (y ,n) gets transformed to U.' 0i (y ,ft) 

I40/0, ft) = Ui(y , n)e iaab ^ n -yM(yo, ft + i)U}(y + 1, n)Ul (Vo, ft) (43) 
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ab ^- y M(yo,n + l)Uj(y + l,n) 



or 



Ui(y + 1, ft) - Ui(y + l,fi)U&(y , n + i^"^"^^, n + i) (44) 



which is a multiplication of the link variable by an SU (N) matrix and it can 
be reabsorbed in a change of variables. However C/j(yo + 1,^) also appears in 
the plaquettes 11^- (y + 1,^) and n 0i (y + l,n). Up to terms 0(a 2 ), the net 
effect will be that in Hij(yo + 1, ft), 



i.e. that a monopole has been added at y + 1, and that a change like the one 
in Eq. (43) is produced at time y + 1. 

By successive iteration of this procedure one eventually come at a time y' 
where an anti-monopole of type a is situated, and then the procedure stops. 

At T = the correlator (p,(x,t)/i(x,0)) can be measured, and by cluster 



whence | (/x) | can be extracted. 

At T 7^ ([/,) is measured directly, and C* boundary conditions in time are 
needed [1,3,4]. 

By the appropriate choice of bf (x — y) in Eq. (35) a generic number of 
monopoles and anti-monopoles can be created at time yo- 

In numerical simulations it is convenient to measure 




property 



(fl(x,t)f,(x,0)) ~ |(/i)| 2 + ce 

t— >oo 



Mt 



d 



log(/i) 




which is much less noisy, and in terms of which 
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The statement that (p) 7^ for T < T c in the infinite volume limit corresponds 
to have p volume independent and finite at large volumes, (p) = 0, for T > T c , 
is obtained if 

p -> -JfeL + fc' (Jfe > 0) (45) 

as the spatial size of the lattice, L, diverges. A behavior like Eq. (45) is numer- 
ically easy to test, and means that {p} is strictly zero in the thermodynamical 
limit. A direct measurement of (p) would only produce a value which is zero 
within (large) errors. 



4 Conclusions 

We have analyzed how to investigate dual superconductivity of the vacuum in 
the confined phase of SU(N) gauge theories for arbitrary N, by deriving in 
detail the abelian projection, its symmetry properties, and the construction 
of a disorder parameter. Numerical simulations are in progress. We plan to 
measure (p), or better p, as a function of /3 — (3/N 2 for different values of 
N. If the ideas about 1/N apply pN 2 should be a universal function of (3, at 
sufficiently large N. Also the comparison of the different choices of (p) and 
of different abelian projections is of interest, on the way to understand the 
nature of dual excitations of QCD. 

This work is partially supported by MIUR, project "Teoria e Fenomenologia 
delle particelle elementari", and by INTAS, Project 00-0110. B.L. is supported 
by the Marie Curie Fellowship No. HPMF-CT-2001-01131. LDD thanks the 
Physics Department in Oxford for hospitality during the final stage of this 
work. 



5 Appendix 

We want to show that the general diagonal, x independent solution of Eq. (1) 
has the form of Eq. (9). Let us call Hi (i = 1. . . N — 1) the independent 
generators belonging to the Cartan algebra of the group, E a the generators 
belonging to the root a. 

The Lie algebra reads then 

[H i: 11^=0 [H i ,E a ] = a i E a (46) 
[E a , Ep] = N a ^pE a+ p [E a , E_ a ] = otiHi 
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In the fundamental representation 



El = E„ a Tr{HiE a } = 
In general for a diagonal $ 

A, = A* Hi + A«E a 



(47) 
(48) 



We want to to solve Eq. (9) 

X 2 (A^, 0) = Tr (0L4 M , A v \) + Tr (0 [L4 M , 0], [A,, 0]]) = (49) 

Using the multiplication law of the algebra, Eq. (46), and Eq. (48) 

[A„ A u ] = A^A? [E a , E p \ + A^ ai E a - A^ ai E a (50) 

The only contribution to X 2 comes from the first term of Eq. (49) and gives 

-r rv-A a A~ a 

As for the second term of Eq. (42), since 



c i AZ[H i ,E a ]=c i a i A°E a 



one has: 



[$, AJ, [$, = A^cmcjPjiE^ E p 



and the second contribution to X 2 is 

1 



A'^A l , a (ciai) 3 



In summary X 2 = is equivalent to 



Va 



(51) 



13 



By gauge transformations we can order the eigenvalues of $ in decreasing 
order 



$ = diag(y?i . . . ip n ) ip ± > ip 2 > . . . <p r , 



(52) 



It is easy to find N — 1 independent solutions of Eq. (51), & g = 2fi g where fjfl 
are the fundamental weights 



p p p q q q 

N' N' ' ' ' ' N' ~N' ~N' ' ' ' ' ~N 



V 



/ 



(53) 



with q = 1,2,..., N- 1. 
Indeed for the simple roots 



2Tr {/j a a b } = 5 ab 



A generic positive root a is the sum of simple roots, so that 

2Tr{/x a o;} = 1 if a contains a a 
2Tr{/i a a} = if a does not contains a a 

The sign changes for negative roots. The /z's form a complete set for diagonal 
operators. The general solution will be of the form 

$ = 2^cV 

A 

2Tr($a A ) = c A = 0, ±1 
it is easy to see from Eq. (53) that 



^ ^ c i 
1 2 N 

* * C 2 
$2 - $ 3 = ^ 



Since the eigenvalues are ordered all the coefficients c A must be non negative, 
if two c's were different from zero, c l , c? , then for the root a = a 1 + a l+1 + . . . a? 
we would have for the scalar product 

2Tr = 2 
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which is impossible. Therefore only one c can be different from zero and $ a 
provides the generic solution of Eq. (9). 
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